TECHNICAL 

LIBRARY 

|ad.<^/3I  /5) 


TECHNICAL  REPORT  ARBRL-TR-02499 


TRANSIENT  COMBUSTION: 
MODELS  VS  ONERA  DATA 

Carl  W.  Nelson 

June  1983 


US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 
BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND,  MARYLAND 


Approved  for  public  release;  distribution  unlimited. 


DnC  QUAUT?  INSPF.tTrHSB  ft 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Additional  copies  o£  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.  S.  Department  of  Commerce,  Springfield,  Virginia 

22161. 


The  findings  in  this  report  are  not  to  be  construed  as 
an  official  Department  of  the  Army  position,  unless 
so  designated  by  other  authorized  documents. 


The  use  of  trade  names  or  manufaoturera '  names  in  this  report 
does  not  aonatitute  indorsement  of  any  aormeroial  produot. 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

TECHNICAL  REPORT  ARBRL-TR-02499 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT’S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitie) 

TRANSIENT  COMBUSTION:  MODELS  VS  ONERA  DATA 

S.  TYPE  OF  REPORT  8  PERIOD  COVERED 

Final 

6.  performing  ORG.  report  NUMBER 

7.  AUTHORf#;  j 

Carl  W.  Nelson 

a.  CONTRACT  OR  GRANT  NUMBER('eJ 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Ballistic  Research  laboratory 

ATTN;  DRDAR-BLI 

Aberdeen  Proving  Ground,  MD  21005 

10.  PROGRAM  ELEMENT,  PROJECT.  TASK 

AREA  &  WORK  UNIT  NUMBERS 

RDT&E  1L161102AH43 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Armament  Research  &  Development  Command 

US  Army  Ballistic  Research  Laboratory  (DRDAR-BLA-S) 
Aberdeen  Proving  Ground,  MD  21005 

12.  REPORT  DATE 

June  1983 

13.  NUMBER  OF  PAGES 

31 

14.  MONITORING  AGENCY  NAME  &  ADDRESS(if  different  from  Controiiing  Office) 

IS.  SECURITY  Cl.ASS.  (ol  thi*  report) 

UNCLASSIFIED 

ISa.  declassificAtion/downgrading 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  thie  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abetract  entered  in  Bfock  20,  if  different  fro3i  Report) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  retrera*  aide  if  neceaaaiy  and  identify  by  6/oclr  number) 

Propellant  Computer 

Combustion  Transient  Combustion 

Internal  Ballistics 
Modeling 

20.  ABSTRACT  (Cdatttxum  aa  reverma  aldk  ff  nacaaemry  and  ideniify  by  6/oclr  number) 

meg 

Models  of  transient  combustion  of  solid  propellants  were  used  to  try  to 
reproduce  transient  regression  rates  measured  in  experiments  at  ONERA 
(French  space  agency).  Variations  of  the  thermal  theory  approach  could  not 
reproduce  the  measured  rates  with  "reasonable"  assignment  of  unmeasurable 
parameters  in  the  models.  Attempts  were  made  with  the  Kooker-Zinn, 
Zeldovich,  and  flame  sheet  approaches  to  the  gas  phase  heat  feedback. 
Calculated  transients  were  smaller  than  the  measurements.  (Cont’d) 


EX) 


FORM 
I  JAM  73 


1473 


EDrUOM  OF  1  NOV  65  IS  OBSOLETE 


UNCLASSXEIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whert  Data  Entered) 


UNCLASSIFIED 


TABLE  OF  CONTENTS 

Page 

List  of  Figures . ...5 

INTRODUCTION . 7 

THE  MODELS . 8 

NUMERICAL  METHODS . . 

THE  ONERA  EXPERIMENT . 11 

RESULTS . 14 

DISCUSSION . 18 

REFERENCES . ,21 

LIST  OF  SYMBOLS . 23 

DISTRIBUTION  LIST . 25 


LIST  OF  FIGURES 


Figure  Page 

1  ONFRA  Apparatus . 12 

2  Regression  Rate,  Test  7.. . «.12 

3  Regression  Rate,  Test  9 . 12 

4  Relative  Regression  Rates ,  ONERA . . . 1^ 

5  Effect  of  Surface  Heat  Release,  Kooker  Model . 15 

6  Effect  of  Activation  Energy,  Kooker  Model . 15 

7  Transition  Boundary,  Kooker  Model . 15 

8  Unstable  Case,  Constant  Pressure . 16 

9  Effect  of  Temperature  Sensitivity,  Zeldovich . 17 

10  Effect  of  Activation  Energy,  Zeldovich,; . .17 

11  Zeldovich  Stability  Diagram . 19 


5 


!•  INTRODUCTION 


Solid  propellant  burning  rate  is  usually  assumed  to  depend  only  on 
pressure  for  the  convenience  of  the  design  engineer*  There  remains,  however, 
the  naggging  suspicion  that  things  are  probably  not  that  simple*  Although 
fast  digital  computers  have  obviated  the  need  for  simple  algebraic  expressions 
(r=ap^),  we  have  given  only  a  little  attention  to  the  possibility  that  the 
burning  rate  dependence  may  extend  to  other  variables  such  as  rate  of 
pressurization,  pressure  history,  or  stress  history  in  the  solid*  In  high 
performance  applications,  burning  rate  may  then  be  fully  transient  (read: 
quasi-steady  treatment  is  inadequate)* 

Even  if  transient  combustion  exists,  however,  its  effect  on  internal 
ballistics  is  still  in  doubt*  Catastrophies  of  overpressure  have  had  at  least 
one  other  credible  explanation:  igniter  operation,  slow  ignition,  bed 
permeability,  grain  fracture*  Any  prediction  of  transient  burning  has  had  a 
serious  limitation:  quasi-steady  flame  in  rapid  pressure  change,  arbitrarily 
simple  (or  ignored)  chemistry*  Experimental  evidence  of  transient  burning  is 
practically  nonexistent.  And  without  demonstrable  evidence  of  a  measurable 
effect,  the  theories  will  not  compete  with  the  ignition  and  mechanical 
explanations • 

Predictions  of  transient  combustion  have  come  from  models*  One 
experiment  with  clear  results  was  conducted  by  the  French  ONERA  (Office 
National  d* Etudes  et  de  Recherches  Aerospatiales),^  The  pressure  changes  were 
easily  measured,  the  data  reduction  was  uncomplicated,  propellant  properties 
were  reported,  and  the  results  were  reproduced  in  tests  with  more  than  one 
propellant*  It  is  an  experiment  against  which  models  can  be  measured* 

Models  are  many*  Unfortunately,  the  few  that  address  the  monotonic 
sharply  rising  pressure  do  not  agree*  The  most  popular  approach  is  the 
thermal  theory  where  the  thermal  condition  of  the  solid  dominates  the  burning 
rate  response  which  depends  on  surface  tempeT^ature * 

2  o 

Krier  pioneered  the  thermal  theory  approach*  Kooker  and  Nelson"^  used  it 
to  calculate  sharp  transient  responses  in  a  gun  pressurization  and  Kooker^ 
used  it  for  unsteady  rocket  motor  calculations.  All  assumed  uniform  spatial 
heat  release  to  integrate  the  flame  energy  equation  to  determine  the  heat 


^e7.  Bvulard,  P*  Kuentzmanriy  and  R.  Kling,  ^^Reponse  d^un  Propevgol  Solide  a  un 
Echelon  de  Pression,”  La  Recherche  Aerospatiale,  Vol.  5,  pp.  279-^287,  1975. 

2  , 

H.  Kr^er,  J.  T^ien,  W.  Sirignano,  and  M.  Surmerfield,  "Nonsteady  Burning 
Phenomena  of  Solid  Propellants:  Theory  and  Experiments,"  AIAAJ,  Vol.  6. 
pp.  278-285,  1968. 

D.  Kookev  and  C.  Nelson,  "Numevieal  Solution  of  Solid  Propellant  Transient 
Combustion, "  ASMS  J  Heat  Transfer.  Vol.  101,  pp.  359-364,  1979. 

4 

D.  Kooker  and  B.  Zinn,  "Triggering  Axial  Instabilities  in  Solid  Rockets: 
Numerical  Predictions,  "  AIAA  Paper  73-1298,  1973. 
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feedback  at  the  propellant  surface.  Sununerfield ,  et  al,^  used  the  Zeldovlch 
heat  feedback  approach  which  avoided  any  specification  of  the  flame 
structure.  Nelson^  showed  that  for  the  same  gun  pressurization  rate,  the 
Krier  approach  and  the  Zeldovich  approach  predicted  quite  different 
transients. 

Kumar  and  Kuo^  used  the  Zeldovich  approach  to  predict  a  runaway 
regression  rate  for  constant  pressurization  rate.  In  contrast,  Nelson”  and 
Caveny®  found  only  modest  excursions. 

This  report  will  apply  such  models  to  the  ONERA  transient  burning 
experiment . 


II.  THE  MODELS 

Thermal  theory  models  all  solve  the  transient  heat  conduction  problem  in 
a  semi-infinite,  inert  solid 

r^T  .  8T  ,  a  3T  .  .  /  /  n  /IN 

^  J  ^  ^  °  ’  -<x<0  .(1) 

In  the  convective  terra,  the  regression  rate  (r)  depends  on  the  surface 
temperature  (Tg)  which  makes  the  equation  nonlinear  and  analytically 
intractable. 


Usually  the  thermal  properties  are  assumed  constant  and  the  equation 
reduces  to 


aT 

at 


+ 


(2) 


Allowing  temperature  dependence  of  the  properties  produces  no  startling 
deviations; 


W.  Summer  fields  L.  Caveny^  /?•  Battista^  N*  Kubota^  Y*  Gostinsev  ^  and 
H.  leoda^  ^Theory  of  Dynamic  Extinguishment  of  Solid  Propellants  Dith  Special 
Reference  to  Ronsteady  Heat  Feedback  Law/*  J  Spacecraft  and  Rockets,  Vol,  8, 
pp.  251^258,  1971. 

C,  nelson,  "Response  of  Three  Types  of  Transient  Combustion  Models  to  Gun 
Pressurization, "  Combustion  and  Flame,  Vol.  32,  pp.  317-319,  1978. 

n 

M,  Kumar  and  K.  Kuo,  "Dynamic  Burning  Effects  in  the  Combustion  of  Solid 
Propellants  with  Cracks  and  the  Use  of  Granular  Bed  Combustion  Models, "  Naval 
Weapons  Center  Report  TP-6193,  1980. 

8 

L.  Caveny,  M.  Suimerfield,  and  C.  Nelson,  "Ignition  Transients  and 
Pressurization  in  Closed  Chambers,"  BRL  Memorandum  Import  2558,  1975  (AD  A017747). 

^C.  Nelson,  "Transient  Combustion  Calculations  with  Variable  Thermal 
Properties,"  BRL  Report  ARBRL-TR-02294,  1981  (AD  A098657). 
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The  regression  rate  is  assumed  to  depend  only  on  surface  temperature  in 
an  exponential  with  two  constants 

r  -  «p  (  -  E^/R„  T^)  .  (3) 

This  commonly  assumed  form  fits  either  a  sublimation  or  a  pyrolytic 
gasification  but  details  of  the  gasification  are  too  poorly  known  to  assign 
firm  values  to  the  two  constants,  pre-exponential  (A^)  and  activation  energy 
(Eg).  Surface  temperature  varies  continuously  and  comes  from  the  solution  of 
the  conduction  equation. 


The  initial  condition  is  a  given  regression  rate  and  a  given 

temperature  distribution  in  the  solid.  A  frequent  condition  is  steady  state 
burning  with  constant  propellant  properties  for  which  the  initial  condition  is 
the  well  known 


T  =  Tq  +  (Tg-Tg)  exp  (rx/oc) 


(4) 


Other  conditions  occur  occasionally  such  as  variable  properties^  or  cold 


solid 


The  boundary  condition  far  into  the  semi-infinite  solid  is  the  constant 
initial  temperature , 


TC-^,  t)  =  T 


(5) 


The  surface  boundary  condition  distinguishes  the  models.  Heat  feedback 
from  the  flame  and  from  any  exothermic  reaction  at  the  solid  surface  combine 
to  form  a  gradient  boundary  condition 


3  T 

^Ig  =  g  (  r,  Tg,  p,  ...)  . 


(6) 


The  Krier  type  condition  of  uniform  heat  release  in  the  flame  is  typified 
by  the  Kooker  model^  where 


ill 

3x '  s 


+  Z^Q^p^/r 


(7) 


and  Z2  are  constants  which  depend  on  the  thermal  and  physical  properties, 


Q 

The  Zeldovich  approach^  leads  to 

|1|  =  -  [T  -  T  -  -  In  (r/r  )] 

dx's  a  s  o  a  s 


(8) 


A  flame  sheet  assumption  where  all  the  heat  release  occurs  at  the  final 
flame  temperature  is  common  in  nitrocellulose  base  propellants  studies.  The 
boundary  condition  there  is 


ill 

3x '  s 


m 


=  Z3r  (Tg  -  T^)  +  Z^  p  /r 


(9) 


Z^  and  Z^  are  also  constants  depending  on  physical  properties. 


9 


Note  that  both  Eq.  (7)  and  Eq.  (9)  assume  that  the  heat  release  term  In 
the  gas  phase  depends  on  the  Instantaneous  pressure 

w  =  Wq  p®  .  (10) 

The  cause  for  change  Is  a  changing  pressure  which  is  Imposed  on  the 
combustion.  In  some  applications  the  pressure  may  come  from  a  coupled 
solution  of  the  chamber  conditions.  In  the  present  application  it  is  simply 
Imposed  as  the  ONERA  measured  pressure  history. 


III.  NUMERICAL  METHODS 

It  is  first  convenient  to  transform  the  problem  into  nondimensional 
variables.  Kooker's  transformation^  Is  convenient. 


T  =  t  (r^/a)  . 

(11) 

n  =  X  (r  /a)  . 
o 

(12) 

0  =  (T-T)/(T  -T). 

o  so  o 

(13) 

H  =  Q  /  [c  (T  -  T  )]  . 
s  s  so  o 

(14) 

E  =  E  /(R  T  )  . 

8  U  SO 

(15) 

A  =  E  (1  -  T  /T  )  . 

o  so 

(16) 

p  =  p/p_^ . 

(17) 

The  new  equations  are 

30  30  32© 

IT  °  O.and 

(18) 

R  =  exp  [A(0  -1)  /  (1  +  A/E  (0  -1))  . 

s  s 

(19) 

The 

initial  condition  is 

0(n  ,0)  =  exp  (Rn ) 

(20) 

The 

cold  boundary  condition  Is 

0  (-«,  T)  =  0  . 

(21) 

The  surface  boundary  condition  for  each  model  Is 

as  follows : 

1 .  Kooker 

30 

R[H  +  (0^-1)  (1-  Cp/Cg)]  +  (1-H)  P/r"*  . 

(22) 

10 


2. 


Zeldovlch 


^  ■  «  ‘®s  - 

p  SO  o 

3.  Flame  Sheet 

=  R  [H  +  (0-1)  -  (c  /c  )(0.-l)]  +  Z,P“/R  .  (24) 

oTl  S  p  8  I  J 

Is  another  constant  depending  on  physical  properties • 

The  numerical  solution  used  an  Implicit  method  with  a  tridiagonal  matrix 
Inversion.  Kooker  and  Nelson^  showed  this  method  was  competitive  with  other 
Implicit  methods  although  It  suffers  a  small  error  when  the  grid  spacing 
varies  greatly  from  one  end  of  the  grid  to  another,  but  the  resultant  accuracy 
loss  Is  small. 

Numerical  accuracy  was  aided  by  using  a  three-point  backward  difference 
scheme  at  the  surface.  Accuracy  and  efficiency  were  combined  by  using  a  grid 
spacing  which  varied  Inversely  with  the  expected  temperature  gradient,  small 
at  the  surface  and  growing  with  distance  Into  the  Interior.  The  surface  grid 
spacing  was  varied  over  the  range  .0005<n<. 023  with  no  useful  effect  on 
accuracy.  Errors  of  10  to  10  In  the  residual  of  Equation  (18)  were  small 
compared  to  the  smallest  term  In  the  conduction  equation  which  was  of  order 
one.  With  typical  values  the  smallest  grid  spacing  allowed  cell-to-cell 
temperature  differences  of  only  about  0.1  K  at  the  surface.  The  combination 
of  time  step  and  grid  spacing  produced  values  of  the  stability  parameter 
(At/CAti^))  ranging  from  0.05  to  200  with  little  effect.  (An  explicit  scheme 
must  not  exceed  0.5.) 

Equation  (18)  was  linearized  by  guessing  the  regression  rate  and 
Iterating  until  successive  solutions  converged  In  both  surface  temperature 
(0^)  and  regression  rate  (R).  Convergence  criteria  were  lO"^  for  the  surface 
temperature  and  lO”  for  the  regression  rate  (both  of  order  one).  Relaxing 
that  convergence  criterion  to  as  large  as  lO”^  had  no  appreciable  effect  on 
the  solution  when  results  were  compared  with  ONERA  combustion. 


IV.  THE  ONERA  EXPERIMENT 

The  ONERA  apparatus  Is  shown  In  Figure  1.  Two  Independent  steady 
pressure  combustion  chambers  are  separated  by  a  passage  blocked  by  a  shear 
disc.  In  Figure  1  the  white  areas  are  the  combustion  gases,  and  the  speckled 
areas  are  the  burning  solid  propellant.  The  propellant  geometry  Is 
unimportant.  When  both  chambers  reach  steady  pressure,  the  shear  disc  Is 
ruptured  to  produce  two  experiments  -  a  pressurization  and  a 
depressurization.  The  depressurization  chamber  starts  at  about  70  bars,  and 
the  pressurization  chamber  at  about  20  bars.  (Figures  1  through  4  are 
reprinted  from  the  ONERA  paper, ^) 


\ 

/ 

) 

i 
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70  bars 


bars 


Figure  1,  ONERA  Apparatus 


Pressure  change  rates  were  low  by  gun  standards, 
pressure  chamber  was  about  MPa/s  compared  to  about 
gun.^ 


The  rate  in  the  low 
10°  MPa/s  in  a  tank 


A  quasi-steady  analysis  reduced  the  pressure  time  history  to 
instantaneous  regression  rate.  Typical  plots  of  regression  rate-time  is  shown 
in  Figures  2  and  3  for  a  positive  pressurization.  The  line  v,  shows  the 
quasi-steady  pressure  dependent  regression  rate. 


Figure  2.  Regression  Rate,  Test  7 


Figure  3.  Regression  Rate,  Test  9 
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Data  after  about  3  ms  should  be  ignored  because  the  assumed  choked  flow 
through  the  passage  no  longer  applies. 

The  figure  of  merit  in  the  experiment  is  relative  regression  rate 
[instantaneous  rate  divided  by  quasi-steady  rate  (v,)].  Figure  4  shows  the 
time  dependence  of  the  relative  rate  in  the  pressurization  chamber  for  three 
distinct  tests.  It  is  this  result  that  the  models  seek  to  duplicate.  Note 
that  the  dotted  curve  on  Figure  4  is  not  an  experimental  result;  it  should  be 
ignored  in  the  context  of  this  paper.  Note  further  that  there  is  some 
discrepancy  between  the  plots  in  Figures  2  and  3  and  in  Figure  4.  Figure  3 
particularly  is  not  typified  by  the  curves  of  Figure  4.  Nevertheless,  Figure 
4  is  the  reported  result. 


The  propellants  were  composites, 
as  reported  by  Brulard^  are: 


Their  physical  and  chemical  properties 


Density,  kg/m^,  p  1735 

Heat  Capacity,  J/kg  K,  Cg  1318 

Initial  Temperature,  K,  T^  293  ^ 

Conductivity,  J/mKs,  X  0.46 

Temperature  sensitivity,  %/K,  0  0.15 

Surface  activation  energy,  kcalFmol,  Eg  25.83 

Reference  surface  temperature,  K,  Tg^  1000 

Regression  rate  pressure  exponent,  n  0.15 


The  physical  and  thermal  properties  can  easily  be  accepted.  The 
combustion  properties  are  open  to  more  debate.  Temperature  sensitivity  and 
pressure  exponent  can  be  directly  measured  and  thus  accepted.  Activation 
energy  and  surface  temperature  are  subject  to  doubt.  Direct  measurement  of 
the  activation  energy  would  require  measurement  of  a  12K  surface  temperature 
change  as  pressure  increased  from  25  to  80  bars.  (Regression  rate  increases 
from  6  to  7  mm/s.)  Whether  the  reference  surface  temperature  is  accurate  is 
of  only  minor  importance  for  the  models;  the  changes  from  it  dominate 
regression  rate  transients  at  such  a  high  activation  energy. 

Some  models  also  need  the  surface  heat  release  to  obtain  the  non- 
dimensional  variable  (H).  For  double  base  propellants  Kubota  inferred 
values  of  0.7  to  0.8  for  H  varying  inversely  with  regression  rate.  Ko^ker^ 
used  values  in  the  same  range  based  on  inference  from  Becks tead,  et  al. 


*  Private  commimieation,  J.  Brulard,  1976. 

Kubota,  T.  Ohlmillev,  L.  Caveny,  and  M.  Sumnerfield,  ^^Site  and  Mode  of 
Action  of  Platonization  in  Double  Base  Propellants/^  AIAAJ,  Vol.  12, 
pp.  1709-1714,  1974. 

Becketead,  R.  Derr,  and  C.  Price, **A  Model  of  Composite  Solid  Propellant 
Combustion  Based  on  Multiple  Flames,^*  AlAAsL  Vol.  8,  pp.  2200-2207,  1970. 
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Fi{.  1  3.  Comparaison  de  diverses  repcnses  k  I’^chelon  posicif. 

- Combustion  radiale  (composition  sans  aluminium). 

-  Combustion  frontale  1  et  2. 

— R^f^rence  [HJ- 


Figure  4.  Relative  Regression  Rates,  ONERA 


V.  RESULTS 

The  models  do  not  reproduce  the  ONERA  regression  rates.  No  calculated 
excursion  even  doubled  the  quasi-steady  regression  rate  which  the  experiment 
found  to  be  3-4  times  the  quasi-steady  rate. 

Kooker  Model 


The  Kooker  model  predicted  only  a  50%  excursion  for  the  nominal 
propellant  properties.  Higher  activation  energies  or  higher  surface  heat 
release  typically  enlivens  responses.  Pushing  both  higher  did,  indeed, 
produce  a  sharp  excursion  but  of  questionable  credibility.  The  peak 
regression  rate  was  either  less  than  2  or  more  than  30.  The  jump  occurred  in 
a  narrow  range  with  no  apparent  middle  ground.  Figures  5  and  6  show  the 
jump.  Note  the  ordinate  is  nondimensional  rate  (instant  rate/quasi-steady 
rate).  ^ 


7  P 

D.  Kooker  and  C.  Kelson,  "Numeriaal  Solution  of  Three  Solid  Propellant 
Combustion  Models  During  a  Gun  Pressure  Transient,"  BEL  Report  1952,  1977 
(AD  A025250). 
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KOOKER  MODEL 

0.75 


Figure  5.  Effect  of  Surface  Heat 
Release ,  Kooker  Model 


KOOKER  MODEL 


Figure  6*  Effect  of  Activation 
Energy ,  Kooker  Model 


Figure  7  shows  the  combination  of  Eg  and  H  which  forms  the  boundary 
between  high  and  low  excursion.  The  sharp  transition,  if  physical,  would  bode 
ill  because  reasonable  experiments  cannot  distinguish  such  fine  differences  in 
either  variable. 


This  behavior  contrasts  with  earlier  findings that 
rate  varies  continuously  with  heat  release  although  those 
lower  activation  energy  (15  kcal/mol)  and  faster  pressure 


the  peak  relative 
findings  were  at  a 
rise  (10^  MPa/s) . 


Figure  7.  Transition  Boundary,  Kooker  Model 
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An  explanation  arises  from  linear  stability  consideration*  Kooker 
showed  that  some  initial  conditions  are  inherently  unstable*  With  no 
disturbance  to  the  initial  condition  (like  pressure  change),  a  spiking 
response  can^be  generated  merely  by  a  disturbance  as  small  as  errors  in  the 
computation*  And,  indeed,  exactly  that  happened;  see  Figure  8  where  constant 
pressure  produced  spikes  for  nominal  (25830)  and  high  H  (0*8)  -  a  point  in 
the  high  response  region  of  Figure  7* 

Figure  8  says  that  calculations  of  transients  in  the  high  response  region 
of  Figure  7  cannot  be  trusted* 

Zeldovich  Model 


The  Zeldovich  approach  calculated  no  transient  over  about  20% * 

Increasing  the  activation  energy  (Eg)  to  as  much  as  50  kcal/mol  could  not 
induce  a  great  excursion,  either* 

The  parameter  which  most  rapidly  affects  transients  in  the  Zeldovich 
model  is  the  temperature  sensitivity*  The  heat  feedback  has  an  excursion 
correction  term  which  is  inversely  proportional  to  the  temperature 
sensitivity.  If  temperature  sensitivity  is  treated  parametrically  to  test  the 
effect  of  its  variation.  Figure  9  results. 


^Thanks  to  D.E.  Kooker  for  pointing  out  the  poeeibility  of  the  linear 
inetahility. 


16 


ZN  MODEL 


TIME  (ms) 

Figure  9.  Effect  of  Temperature  Sensitivity,  Zeldovich 


At  high  temperature  sensitivity  the  excursions  can  be  created;  but  they 
are  followed  by  extinguishment,  which  the  propellant  did  not  see. 

The  same  behavior  occurs  when  the  temperature  sensitivity  is  raised,  and 
activation  energy  is  varied.  Figure  10  shows  that  result. 


ZELDOVICH  MODEL 


TIME  (ms) 

Figure  10.  Effect  of  Activation  Energy,  Zeldovich 
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Flame  Sheet  Model 


The  flame  sheet  model  calculated  no  excursions  over  2%.  Such  a  result 
could  have  been  foreseen  by  noting  that  when  the  flame  reaction  rate  varies 
with  the  two  terms  in  the  heat  feedback  boundary  condition  of  Equation  (9) 
cancel,  and  steady  state  is  assured. 


VI,  DISCUSSION 

The  models  failed  to  reproduce  the  experimental  results.  Which  is  wrong 
-  model  or  experiment?  Being  unwilling  to  challenge  the  data,  we  will  assume 
something  is  amiss  with  the  models, 

Quasi-Steady  Flame 

The  concept  of  an  inert  solid  allows  the  heat  conduction  to  govern  the 
solid  interior  temperature  profile.  Propellants  with  high  enough  activation 
energy  to  be  stable  before  ignition  would  be  inert  except  close  to  the  surface 
if  heat  conduction  is  the  only  mechanism  of  energy  transfer  to  the  unreacted 
solid.  The  most  likely  candidate  for  error,  then,  is  the  surface  heat  release 
or  the  flame. 


The  quasi-steady  flame  is  only  a  convenience.  It  depends  on  an  argument 
that  the  characteristic  time  for  pressure  change  (,001  sec  in  this  experiment) 
is  larger  than  the  flame  characteristic  time  (order  10"^  sec) ,  as  in 
Summerf ield,^  In  a  gun  problem  where  the  pressure  characteristic  time  is  , 

closer  to  10""^  sec,  the  quasi-steady  flame  is  less  appealing  although  it  is 
still  fast  with  respect  to  the  solid.  But  for  the  ONERA  experiment  it  seems  a 
justifiable  assumption. 

The  form  of  the  steady  gas  phase  chemistry  may  be  incorrect.  Neither  of 
the  two  extremes,  uniform  or  flame  sheet,  gave  the  correct  answer.  Perhaps 
some  intermediate  approach,  or  some  combination  approach  like  Merkle,  et  al, 
could  help.  At  least  it  adds  more  adjustable  constants,  Miller^s  idea^  of  a 
quasi-constant  heat  release  is  an  alternate  version  of  a  uniform  release. 


Using  an  unsteady  gas  phase  would  add  complications,  A  description  of 
the  flame  would  require  kinetics  and  transport  information  not  generally 
available  and  a  substantial  increase  in  computing,  Suhas  and  Bose*^  attempted 
a  less  rigorous  step  by  postulating  that  transients  in  blowing  increased  the 


Merkle,  5.  Turk,  and  M.  Sutmerfield,  "Extinguishment  of  Solid  Propellants 
by  Depressurization:  Effects  of  Propellant  Parameters,"  AIAA  Paper  69-176, 
1969, 

Miller,  "An  Idealized  Model  of  Homogeneous  Solid  Propellant  Combustion," 
Combustion  and  Flame.  Vol,  46,  pp,  51-73,  1982. 

1 5 

H.  Suhas  and  T.  Bose,  "A  Mathematical  Model  to  Predict  Transient  Burning 
Rate  and  Decay  Rates  for  Extinction  of  Composite  Propellants , "  Combustion 
and  Flame.  Vol*  28,  pp.  145-153,  1977* 
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heat  transfer  to  the  surface  in  direct  proportion  to  the  blowing;  i,e«,  a 
pressure  decrease  caused  a  regression  rate  increase.  One  consequence  of  that 
approach  was  that  the  positive  pressurization  reduced  the  regression  rate 
which  never  recovered  to  the  initial  rate.  No  positive  excursion  would  result 
from  that  theory. 

Zeldovich 


The  Zeldovich  model  predicts  only  small  excursions  for  this  low 
temperature  sensitivity  propellant.  Earlier  calculations^*  also  produced 
only  low  excursions.  If  the  temperature  sensitivity  could  be  a  variable,  it 
would  be  easier  to  adjust  calculation  to  experiment.  But  it  is  an 
independently  physically  measured  variable  and  not  subject  to  arbitrary 
assignment. 

The  tendency  to  excursions  at  high  temperature  sensitivity  with  Zeldovich 
can  be  seen  in  Zeldovich  stability  diagram  as  shown  by  Stokes.^  For  constant 
Eg,  Figure  11  shows  the  transition  from  stable  to  unstable  burning  as 
temperature  sensitivity  increases  from  0.1%  to  0.4%.  And  this  trend  to 
instability  is  confirmed  in  calculations  of  extinguishment  when  the 
temperature  sensitivity  was  in  the  unstable  region  of  the  diagram. 


K=crp(Ts-Tol 

Figure  11.  Zeldovich  Stability  Diagram 


Only  one  investigator  has  found  large  transients  with  the  Zeldovich 
approach.  Kuo^  calculated  runaway  regression  rates  for  presssurization  at 
10  “10^  atm/s  of  composite  propellant.  Unfortunately,  his  results  could  not 
be  reproduced  with  the  method  described  here  despite  the  fact  that  the 


Stokes^  ^^Application  of  Zeldovich  Heat  Feedback  to  Dynamic  Burning  of 
Throttleable  Fuel^Rich  Propellant  Fuel  Generatore,  Astronauti^ca  Acta, 
Vol.  18,  pp.  395-407,  1973. 
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equations  and  the  implicit  solution  are  the  same.  The  time  and  space 
discretization  are  somewhat  different,  but  variation  over  a  wide  range  of  An 
and  At  produced  no  hint  of  any  substantially  higher  transient.  The  present 
method  can  calculate  spikes  even  as  steep  as  in  Figure  8.  Krier’s  calculated 
runaway"^  was  later  shown  to  be  boundec^  but  it  did  have  a  spike  which  the 
numerical  method  could  calculate.  The  inference  is  that  the  method  would  find 
a  spike  if  it  were  a  solution  to  the  differential  equation  (a  separate 
question  from  whether  it  is  physically  correct). 

An  entirely  different  integration  scheme,  invariant  imbedding,  also 
failed  to  reproduce  Kuo’s  results.  A  second  suspicious  aspect  of  Kuo’s  result 
is  that  the  runaway  occurs  at  about  the  same  pressure  regardless  of  the 
pressurization  rate  (2.4  times  the  initial  pressure).  Unless  Kuo’s  results 
can  be  verified,  it  seems  that  the  Zeldovich  approach  will  not  calculate  high 
transients  for  low  temperature  sensitivity  propellants. 

Initial  Lag 

The  models  do  not  capture  the  initial  drop  in  relative  regression  rate 
shown  in  the  ONERA  results  of  Figure  4.  Such  a  result  is  consistent  with  the 
low  positive  excursions  later  because  other  calculations  with  these  models 
showed  an  initial  relative  regression  rate  lag  as  a  prelude  to  a  positive 
excursion.  Generally,  the  higher  the  positive  excursion,  the  lower  the 
initial  negative  excursion.  Nelson^  found  a  lag  up  to  50%  at  very  high 
pressurization  rate  (7x10^  MPa/s)  with  the  Zeldovich  model  but  only  a  10%  lag 
at  7x10^  MPA/s.  Kooker  found  only  a  small  initial  lag  at  7x10^^  MPa/s.  The 
low  ONERA  pressurization  rate  of  10  MPa/s  should  produce  only  a  small  lag. 

An  initial  lag  is  not  surprising;  the  ideal  thermal  theory  model  assumes  a 
lag.  It  is  possible  that  the  ONERA  data  reduction  exaggerated  the  initial 
response  by  overlooking  some  transient  aspect  of  the  first  compression  wave  to 
enter  the  low  pressure  chamber. 
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ATTN:  B.  Nichols, 

T7,  MS-B284 
P.O.  Box  1663 
Los  Alamos,  NM  87545 
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1  Olin  Corporation 

Smokeless  Powder  Operations 
ATTN :  R.  L.  Cook 
P.O.  Box  222 
St.  Marks,  FL  32355 

1  Paul  Gough  Associates,  Inc. 
ATTN:  P.S.  Gough 

1048  South  Street 
Portsmouth,  NH  03801 

2  Princeton  Combustion 

Research  laboratories 
ATTN:  M.  Summerfield 
N.A.  Messina 

1041  US  Highway  One  North 
Princeton,  NJ  08540 

1  Science  Applications,  Inc. 

ATTN:  H.  S.  Pergament 
1100  State  Road,  Bldg.  N 
Princeton,  NJ  08540 

1  Rockwell  International  Corp. 
Rocketdyne  Division 

ATTN:  J.E.  Flanagan/HB02 
6633  Canoga  Avenue 
Canoga  Park,  CA  91304 

2  Sandia  National  laboratories 
Combustion  Sciences  Dept, 
ATTN:  R.  Cattolica 

D.  Stephenson 
Livermore,  CA  94550 

1  Science  Applications,  Inc, 

ATTN:  R.B,  Edelman 
23146  Cumorah  Crest 
Woodland  Hills,  CA  91364 
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1  Space  Sciences,  Inc. 

ATTN:  M.  Farber 
Monrovia,  CA  91016 

4  SRI  International 

ATTN:  S.  Barker 
D.  Crosley 
D.  Golden 
Tech  Ub 

333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

1  Stevens  Institute  of  Tech 

Davidson  laboratory 
ATTN:  R.  McAlevy,  III 
Hoboken,  NJ  07030 

1  Teledyne  McCorinack-Selph 

ATTN:  C.  leveritt 
3601  Union  Road 
Hollister,  CA  95023 

1  Thiokol  Corporation 

Elkton  Division 
ATTN:  W.N.  Brundige 
P.O.  Box  241 
Elkton,  MD  21921 

3  Thiokol  Corporation 

Huntsville  Division 
ATTN:  D.A.  Flanagan 
Huntsville,  AL  35807 

3  Thiokol  Corporation 

Wasatch  Division 
ATTN:  J.A.  Peterson 
P.O.  Box  524 

Brigham  City,  UT  84302 

1  United  Technologies 
ATTN:  A.C.  Eckbreth 
East  Hartford,  CT  06108 

2  United  Technologies  Corp. 
ATTN:  R.S.  Brown 

R.O.  Mclaren 
P.O.  Box  358 
Sunn5rvale ,  CA  94086 
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1  Universal  Propulsion  Company 

ATTN:  H.J.  McSpadden 
Black  Canyon  Stage  1 
Box  1140 

Phoenix ,  AZ  85029 

1  Veritay  Technology,  Inc. 

ATTN:  E.B.  Fisher 
P.O.  Box  22 

Bowmansville ,  NY  14026 

1  Brigham  Young  University 

Dept,  of  Chemical  Engineering 
ATTN:  M.W.  Beckstead 
Provo,  UT  84601 

I  Director 

Jet  Propulsion  laboratory 

ATTN:  MS  125/159 
4800  Oak  Grove  Drive 
Pasadena,  CA  91105 


1  California  Institute  of 
Technology 

ATTN:  F.E.C.  Culick/ 

MC  301-46 
204  Karman  Lab. 

Pasadena,  CA  91125 

1  University  of  California, 

Berkeley 

Mechanical  Engineering  Dept. 
ATTN;  J.  Daily 
Berkeley,  CA  94720 

1  University  of  California 

Los  Alamos  National  lab. 

ATTN;  T.D.  Butler.  MS  B216 
P.O.  Box  1663 

Los  Alamos,  NM  87545 

2  University  of  California, 

Santa  Barbara 
Quantum  Institute 
ATTN:  K.  Schofield 
M.  Steinberg 

Santa  Barbara,  CA  93106 
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1  University  of  Southern 

California 
Dep t .  of  Cherai s t ry 
ATTN:  S.  Benson 
Los  Angeles,  CA  90007 

1  Case  Western  Reserve  Univ. 

Div.  of  Aerospace  Sciences 
ATTN:  J.  Tien 
Cleveland,  OH  44135 

1  Cornell  University 

Department  of  Chemistry 
ATTN:  E.  Grant 
Baker  Laboratory 
Ithaca,  NY  14853 

1  Univ.  of  Dayton  Rsch  Inst. 

ATTN:  D.  Campbell 
AFRPL/PAP  Stop  24 
Edwards  AFB,  CA  93523 

1  University  of  Florida 
Dept,  of  Chemistry 
ATTN:  J.  Winefordner 
Gainesville,  FL  32601 

3  Georgia  Institute  of 

Technology 
School  of  Aerospace 
Engineering 
ATTN:  E.  Price 
Atlanta,  GA  30332 

2  Georgia  Institute  of 

Technology 
School  of  Aerospace 
Engineering 
ATTN:  W.C.  Strahle 
B.T.  Zinn 

Atlanta,  GA  30332 

1  Hughes  Aircraft  Company 

ATTN:  T.E.  Ward 
8433  Fallbrook  Avenue 
Canoga  Park,  CA  91304 


No.  Of 

Copies  Organization 

1  University  of  Illinois 

Dept,  of  Mech.  Engineering 
ATTN:  H.  Krier 
144  MEB,  1206  W.  Green  St. 
Urbana,  IL  61801 

1  Johns  Hopkins  University/ APL 

Chemical  Propulsion 
Information  Agency 
ATTN:  T.W.  Christian 
Johns  Hopkins  Road 
Laurel,  MD  20707 

1  University  of  Minnesota 

Dept,  of  Mechanical 
Engineering 
ATTN:  E.  Fletcher 
Minneapolis ,  MN  55455 

4  Pennsylvania  State  University 

Applied  Research  Laboratory 
ATTN:  G.M.  Faeth 
K.K.  Kuo 

H.  Palmer 
M.  Micci 

University  Park,  PA  16802 

1  Polytechnic  Institute  of  NY 
ATTN:  S.  lederman 

Route  110 

Farmingdale ,  NY  11735 

2  Princeton  University 
Forres tal  Campus  Library 
ATTN:  K.  Brezinsky 

I.  Glassman 
P.O.  Box  710 
Princeton ,  NJ  08540 

1  Princeton  University 

MAE  Dept.,  Eng.  Quad. 

ATTN:  F.A.  Williams,  D325 
P.O.  Box  710 
Princeton,  NJ  08540 
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2  Purdue  University 
School  of  Aaronautics 

and  Astronautics 
ATTN:  R.  Click 

J.R.  Osborn 
Grissom  Hall 

West  Lafayette,  IN  47907 

3  Purdue  University 
School  of  Mechanical 

Engineering 
ATTN:  N.M.  Laurendeau 
S.N.B.  Murthy 
D.  Sweeney 
TSPC  Chaffee  Hall 
West  Lafayette,  IN  47906 

1  Rensselaer  Polytechnic  Inst. 
Dept,  of  Chemical  Engineering 
ATTN:  A.  Fontijn 

Troy,  NY  12181 

2  Southwest  Research  Institute 
ATTN:  Robert  E.  White 

A.B.  Wenzel 
8500  Culebra  Road 
San  Antonio,  TX  78228 
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1  Stanford  University 
Dept,  of  Mechanical 

Engineering 
ATTN :  R.  Hanson 
Stanford ,  CA  93106 

2  University  of  Texas 
Dept .  of  Chemistry 
ATTN:  W.  Gardiner 

H.  Schaefer 
Austin,  TX  78712 

1  University  of  Utah 

Dept,  of  Chemical  Engineering 

ATTN:  G.  Flandro 

Salt  Lake  City,  UT  84112 

1  Virginia  Polytechnical 

Institute  and 
State  University 
ATTN:  J.A.  Schetz 
Blacksburg ,  VA  24061 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN :  DRXSY-D 

DRXSY-MP,  H.  Cohen 
Cdr,  USATECOM 

ATTN :  DRSTE-TO-F 
Dir,  USACSL,  Bldg.  E3516,  EA 
ATTN:  DRDAR-CLB-PA 
DRDAR-CLN 
DRDAR-CLJ-L 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number^ _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc .) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: _ 

Telephone  Number : _ 


Organization  Address: 


